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ABSTRACT 

Context. Unbiased spectral surveys are powerful tools to study the chemistry and the physics of star forming regions, as they can give 
a complete census of the molecular content and the observed lines probe the physical structure of the source. 

Aims. While unbiased surveys in the millimeter and sub-millimeter ranges observable from ground-based telescopes have previously 
been obtained towards several high mass protostars, very little exists on low mass protostars, believed to resemble our own Sun's 
progenitor. Aiming to fill up this gap, we carried out a complete spectral survey of the bands at 3, 2, 1 and 0.8 mm towards the solar 
type protostar IRAS 16293-2422. 

Methods. The observations covered about 200 GHz and were obtained with the IRAM-30m and JCMT-15 m telescopes during about 
300 hours of observations. Particular attention was devoted to the inter-calibration of the obtained spectra with previous observations. 
All the lines detected with more than 3 cr and free from obvious blending effects were fitted with Gaussians to estimate their basic 
kinematic properties. 

Results. More than 4000 lines were detected (with cr > 3) and identified, yielding a line density of approximatively 20 lines per GHz, 
comparable to previous surveys in massive hot cores. The vast majority (~ 2/3) of the lines are weak and due to complex organic 
molecules. The analysis of the profiles of more than 1000 lines belonging 70 species firmly establishes the presence of two distinct 
velocity components, associated with the two objects, A and B, forming the IRAS 16293-2422 binary system. In the source A, the line 
widths of several species increase with the upper level energy of the transition, a behavior compatible with gas infalling towards a 
~ 1 M object. The source B, which does not show this effect, might have a much lower central mass of ~ 0.1 M G . The difference in 
the rest velocities of both objects is consistent with the hypothesis that the source B rotates around the source A. 
Conclusions. This spectral survey, although obtained with single-dish telescope with a low spatial resolution, allows to separate the 
emission from 2 different components, thanks to the large number of lines detected. The data of the survey are public and can be 
retrieved on the web site http://www-laog.obs.ujf-grenoble.fr/heberges/timasss 
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1. Introduction 

It is well known that the chemical composition of the gas from 
which the star forms influences the process of the star formation 
and it is, in turn, influenced by it. The first obvious example is 



* Based on observations with the IRAM 30 m telescope (IRAM 
is supported by INSU/CNRS (France), MPG (Germany) and IGN 
(Spain)), and with the JCMT 15 m telescope (operated by the Joint 
Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organisation 
of Scientific Research, and the National Research Council of Canada). 



that the Jeans mass depends on the gas temperature which, in 
turn, is set by the molecular line cooling in a large range of den- 
sities and temperatures. The chemical composition of the gas is 
here of paramount importance, as the cooling is indeed domi- 
nated by different species as a function of the gas temperature 
and density and the elemental abundance (Goldsmith 2001 1. A 
second classical example is the slow contraction of the prestellar 
cores which is governed by ambipolar diffusion. Since only ions 
feel the magnetic field which counteracts the gravitational force, 
the chemical composition of the gas, which determines the ions 



2 



Caux et al.: IRAS 16293-2422 spectral survey 



abundance, is crucial. In addition to this, since the gas chemi- 
cal composition is largely affected by the star formation process, 
its study in star forming regions is a powerful diagnostic tool to 
track down the various processes at work. Finally, the study of 
the chemical composition in regions forming solar type stars is 
of particular importance, as it helps to understand the history of 
the formation of our own Solar System. For example, the com- 
parison of the chemical composition of comets with that of solar 
type protostars or protoplanetary disks is used to trace back the 
origin of the former ( |Crovisier et al.|2004 1. The same applies to 
the studies of the molecular content of meteorites. In particular, 
the recent claim that the amino acids found in carbonaceous me- 
teorites may have formed during the first phases of the life of 
the Solar System is based on the measured large deuteration in 
the meteoritic amino acids and in the protostellar environment 
( |Pizzarello & Huang||2005| l. In this context, unbiased spectral 
surveys in millimeter and submillimeter wavelengths are partic- 
ularly relevant as they allow to detect heavy and large molecules, 
and, specifically, complex organic molecules. In summary, un- 
biased spectral surveys in the millimeter to submillimeter wave- 
lengths are a powerful method to characterise the molecular con- 
tent of astrophysical objects, and the only way to obtain a com- 
plete census of the chemical species. 

There are at least two other aspects that make unbiased 
spectral surveys precious tools for studying the star formation 
process. First, in general they provide multiple lines from the 
same molecule, allowing multi-frequency analysis and mod- 
elling. Since different lines from transitions with different upper 
level energies and Einstein coefficients are excited at different 
temperatures and densities, they probe different regions in the 
line of sight. A careful analysis can, therefore, distinguish be- 
tween the various physical components in the beam. If one adds 
also the kinematic information provided by the line profiles, the 
method can be so powerful that it can identify sub-structures in 
the line of sight, even if the spatial resolution of the observations 
is inadequate. In the present article, we provide an example of 
this capability of unbiased spectral surveys. 

Given their powerful diagnostic ability, several unbiased 
spectral surveys in the millimeter and sub-millimeter bands ac- 
cessible from ground have been obtained in the past in the di- 
rection of star forming regions. A complete list of these surveys 
can be found in Herbst & van Dishoeck (2009). By far the most 
targeted sources are hot cores, the regions of high mass proto- 
stars formation, where the dust temperature exceeds the sub- 
limation temperature of the water-ice grain mantles, ~ 100K. 
The combined effect of the mantle sublimation and the high 
gas temperature triggers a singular and rich chemistry. At the 
same time, the relatively high densities (> 107crrT 3 ) and tem- 
peratures (> 100 K) are favorable for the excitation of several 
high lying transitions. The result is that extremely rich spec- 
tra are observed towards the hot cores. About a dozen hot 
cores have been targeted in differ ent bands (|Schilke et al.|1997 



Helmich & van Dis hoeck 1997; Hatc hell et al.||1998t |Schilke 
et al. |2001| |Tercero et al.||2010| i. One of the most studied hot 
cores is the Orion- KL source. Spectral surveys covering almost 
all the bands accessible from the ground have been obtained, 
from about 80 to 900 GHz, detecting thousands of lines from 
hundreds of species and relative isotopologues (see for instance 
Schil ke et al.| <fl997); |Lee et al.| pOUI) ; |Comito et all (|2005); 



Olofsson et al. (2007i and references therein; |Demyk et al 
(|2007|l; ICarvajal et al.[ ( p009] l; |Margules et al.|p009) ; |Tercero 



et al.| ( [20T0| . In addition, t he 500-2000 GHz rang e is observed 



recently launched Herschelp]satellite (Pilbratt et al. 2010|, in the 
Key Program HEXOS (http://www.submm.caltech.edu/hexos/). 
Similarly, the 500 — 2000 GHz range is observed in other hot 
cores, in the Key Program CHESS (http://www-laog.obs.ujf- 
grenoble.fr/lieberges/chess/). Preliminary results of the surveys 
performed in these two Herschel Key Programs can be found in 
(Berg met al.|2010) and (jCeccarelli et a l.|2010| l, respectively. 

Although less massive and less luminous, solar type pro- 
tostars also possess regions where the dust mantles sublimate, 
yielding similar properties as those of hot cores (see |CeccarelIT| 
( |2007| l and references therein). These regions have been baptized 
hot corinos, to make it clear that they share similarities with hot 
cores but are not just scaled version of them (see also |Bottinelli| 
|et al.||2007 1. The interest in observing hot corinos, whose sizes 
are comparable to the Solar System sizes, is amplified by the 
fact that they likely resemble the Solar Nebula. In other words, 
their study corresponds to an archeological study of our ancestor 
Solar System. 

So far, only one (partial) spectral survey has been obtained 
towards a solar type protostar ( Blake et al.|1994 van Dishoeck 



[eFaT1[T995] l. The targeted source was IRAS 16293-24 22 (here 
inafte r IRAS 16293), in the L1689N cloud (d= 120pc, [Loinard| 
|et al.||2"0"08"j ). This survey partially covered the two windows in 
the 200 GHz and 350 GHz bands accessible from ground, and 
has been obtained with the JCMT and CSO telescopes. The 
sensitivity achieved (~ 40 mK ) allowed the detection of 265 
lines from 24 species, namely the most abundant molecules: 
CO, H 2 CO, CH 3 OH, SO, S0 2 etc... Later, more sensitive ob- 
servations have shown that the IRAS16293 line spectrum is rich 
in complex organic mol ecules (jCeccarelli et al.|2000b| |Cazaux 
l et al.|2003) , and doubly ( jCeccarelli et al.|1998| Tand triply ( |Parise 
|et al.|2004 1 deuterated molecules. In retrospective, further sup- 
port for the interest in obtaining an unbiased spectral survey to- 
wards IRAS 16293 is provided by the recent Herschel/HIFI data 
in the 555—636 GHz range: it shows that, while IRAS 16293 has 
much less lines than the 2 x 106 L source NGC6334I, the same 



number of species is detected in both sources (Ceccarelli et al. 
[20T0| . 

Several studies have been carried out towards IRAS 16293, 
both with single dish telescopes and interferometers. The emerg- 
ing overall picture is that IRAS 16293 is a protobinary system 
( |Wootten|1989| |Mundy et al.|1992 1 surrounded by an envelope 
of about 2M Q ( |Crimier et al.||2010| ). The structure of the en- 
velope has been the target of several studies (jCeccarelli et al. 



2000a) jSchoier et al.|2OO2||j0rgensen et al.|20()5| l. The mostre 
cent one by Crimi er et al. ( 2010) > concludes that the envelope 
density follows a r~ z power law at distances larger than about 
1300 AU and r~ 3 ^ 2 innerwards. The grain mantles are predicted 
to sublimate at a distance of 75 AU, where the density is equal to 
2 x 108 crrr 3 . The envelope extends from about 25 AU to about 
6000 AU from the center. Inside the envelope, the two sources, 
A (South-East) and B (North-West), of the binary system are 
separated by about 4" (separation measured from interferome- 
ter observations at a spatial resolution of about 1"), equivalent 
to a linear distance of 480 AU. The source B is brighter than the 
source A in the millimeter continuum and in several "cold enve- 
lope" molecular lines, whereas the source A seems to be brighter 
in several hot corino like molecular fines ( |Kuan et al.||2 004; 
Bottinelli et"aT1[2004} |Chandler et al.||2005) . Finally, |ChandIer] 



et al. ( 2005) > have claimed that source A might be a multiple sys- 



with the HIFI spectrometer ( |de Graauw et al.|2010| l on board the tant participation from NASA 



1 Herschel is an ES A space observatory with science instruments pro- 
vided by European-led principal Investigator consortia and with impor- 
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tem and very recent observations ( Pec h et al.|2010[ ) suggest that 
A is itself a binary system of 0.5 and 1.5 M Q respectively. 

Despite its relatively complex structure at arcsec scales, 
IRAS 16293 remains the brightest and best source to carry out 
a detailed study of the gas chemical composition in the first 
phases of the formation of a solar type star. As discussed above, 
the best way for that is to obtain unbiased spectral surveys, 
as sensitive as possible. In this paper, we present the results 
of the most sensitive unbiased spectral survey of the bands 
between 80 and 366 GHz observable from ground based tele- 
scopes obtained so far in direction of IRAS 16293. This study 
is part of a more general project that plans on observing also 
the 500-2000 GHz frequency range with the spectrometer HIFI 
on board the Herschel satellite, in the context of the Herschel 
Key Program CHESS (Chemical Herschel Surveys of Star form- 
ing regions: http://www-laog.obs.ujf-grenoble.fr/heberges/chess/, 
( |Ceccarellietal.|2010| ). " 

2. Observations 

The observations were obtained at the IRAM-30m (frequency 
range 80-280GHz) and JCMT-15m (frequency range 328- 
366 GHz) telescopes during the period January 2004 to August 
2006. Overall, the observations required a total of about 300 
hours (~ 200 hr at IRAM and ~ 100 hr at JCMT) of observ- 
ing time. The beam of the survey varied between 9" and 33", 
depending on the used telescope and frequency, and the spec- 
tral resolution ranged between 0.3 and 1.25 MHz, correspond- 
ing to velocity resolutions between 0.51 and 2.25 km s~'. The 
achieved rms varied between 4 and 14 mK in 1.5kms _1 bins. 
The observations were centered on the B (North-West) compo- 
nent at a(2000.0) = 16'' 32'" 22.6 m , 5(2000.0) = -24° 28' 33'. A 
and B components, separated by 4", are both inside the beam of 
our observations at all frequencies. However, at the highest fre- 
quencies observed with the IRAM 30 m telescope (i.e. the 1 mm 
band), the attenuation of emission from source A is not negli- 
gible as will be discussed in section 5. All observations were 
performed in Double-Beam-Switch observing mode, with a 90" 
throw. Pointing and focus were checked every two hours on plan- 
ets and on the continuum radio sources 1741-038 and 1730-130. 
TableQ]summarizes the observed bands and the details of the ob- 
servations. More details are given below for the IRAM-30m and 
JCMT- 15 m telescope observations, respectively. Because of the 
different weather conditions during the different runs, the sys- 
tem temperatures varied largely. However, during the data pro- 
cessing, scans with too large system temperatures were removed 
before averaging. 

2.1. IRAM Observations 

The following three bands were almost fully covered by observa- 
tions at the IRAM-30 m telescope: 3 mm band (80-1 15.5 GHz), 
2 mm band (129-177.5 GHz) and 1 mm band (198-281.5 GHz). 

In all IRAM-30 m observations, two frequency ranges were 
observed simultaneously, with two SIS receivers with orthogonal 
polarizations for each frequency range, in the following config- 
uration: 3 mm receivers (A100 & B100) in parallel with 1 mm 
receivers (A230 & B230), and 2mm receivers (C150 & D150) 
in parallel with 1 mm receivers (C270 & D270). Because of the 
limitation (at the time of the observations) of the IRAM-30 m 
backend capabilities in terms of instantaneous frequency band- 
width and spectral resolution, we privileged the largest possible 
spectral bandwidth to cover the IRAM-30 m bands in the small- 
est observing time. For simultaneous observations in the 3 mm 



and 1 mm bands, the VESPA autocorrelator was split in 4 parts, 
2 of them covering the whole IE band of the A&B100 receivers 
(0.5 GHz) with 320 kHz spectral resolution and the two others 
covering half of the IF band of the A&B230 receivers (1 GHz) 
with 1250 kHz spectral resolution. The second half of the IF 
band of the A&B230 receivers was covered with the 1 MHz 
Filter Banks (FB). For simultaneous observations in the 2 mm 
and upper 1 mm bands, the VESPA autocorrelator was split in 4 
parts, 2 of them covering half of the IF band of the C&D150 
receivers (0.5 GHz) with 320 kHz sampling two others cover- 
ing half of the IF band of the C270&D270 receivers (1 GHz) 
with 1250kHz sampling. The second half of the IF band of the 
C&D150 receivers was covered with the 1 MHz FB. 

The configuration for observations in the 2 mm band resulted 
in different spectral resolutions for each half of the IF band of 
the receivers (320kHz ~ 0.65 km s 1 and 1 MHz ~ 2kms~'). 
Therefore, we shifted the tuning frequency of the receivers by 
only 0.5 GHz from one tuning to the next one to cover the en- 
tire 2 mm band at the highest and at the smallest spectral res- 
olution respectively. As a consequence, two different datasets 
were obtained, one at high resolution (generally used for study- 
ing brighter lines), and one at low resolution (for the faint lines). 

2.2. JCMT Observations 

The JCMT-15m observations covered the 328 to 366 GHz fre- 
quency range. They were obtained with 345 GHz SIS receiver 
RxB3 in dual-channel single-sideband (SSB) mode. Each po- 
larization of the receiver was connected to a unit of the ACSIS 
autocorrelator providing a bandwidth of 0.5 GHz for a spectral 
resolution of 625 kHz. At 345 GHz, this yields a velocity resolu- 
tion of about 0.5 km s _1 . 



3. Calibration 

3.1. Method 

At the 30 m telescope, the calibration was performed with a cold 
and a warm absorbers, and the atmospheric opacity was obtained 
using the ATM program (Cernicharo, 1985, IRAM internal re- 
port). At the JCMT, line strengths were calibrated via the chop- 
per wheel method (Kut ner & Ulich|1981) >. 

Our spectral survey does not allow to estimate the calibra- 
tion uncertainties from line observation redundancy: each spec- 
tral range has been observed only once and there is only negligi- 
ble frequency overlap between adjacent spectra. From our simul- 
taneous observations with two receivers in the 1 mm range we 
may estimate the receivers contribution, but to derive the total 
calibration uncertainties of the survey, we have performed a de- 
tailed comparison between our observations and previous obser- 
vations. As our comparisons rely only on observations obtained 
with the same telescopes towards the same position (namely 
source B), there is no bias due to different sources dilution in 
the beams and our results are not affected by the underestimate 
of source A contribution at high frequency. The comparison in- 
cludes virtually all the published data towards IRAS 16293 ob- 
tained with the IRAM-30m and the JCMT-15m telescopes, as 
well as unpublished data obtained with the IRAM-30 m tele- 
scope. The list of the articles used for this comparison is the fol- 
lowing: a) IRAM-30m bands: Ceccarelli et al.|([ l998); Loinard 



etlll ( |2000li ; |Cazaux et al.| p003 ); Wakelam 'etal.|(|2004[ ) ;|Parise 



|van Dishoeck et aL 
( |2002| i; |Parise etal. 



eraL|p002] |2005bb; b) JCMT-15 m band : [Blake et aT](l994|l; 



1995); Loinar cTeTaT1 ( |2000| ); |Schoier et al.| 
2004). Table [Tjreports the percentage of the 
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Table 1. Parameters of the observations at IRAM-30m and JCMT-15m telescopes. 
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IRAM 


80-115.5 


0.32 


0.81-1.17 


VESPA 


2-8 


90-400 


30-21 


0.80-0.78 


11 


28 


15 


IRAM 


129 - 177 


0.32 


0.53-0.72 


VESPA 


5-14 


200-1000 


19-14 


0.76-0.69 


17 


22 


12 


IRAM 




1.0 


1.65-2.25 


1 MHz FB 










10 (6) 


84 (6) 


95(6) 


IRAM 


197 - 265 


1.0 


1.13-1.52 


1 MHz FB 


7-13 


180-1200 


12-9 


0.65-0.51 


17 


36 


10 


IRAM 




1.25 


1.41-1.90 


VESPA 
















IRAM 


265 - 280 


1.25 


1.34-1.41 


VESPA 


9-17 


470-4200 


9 


0.51-0.47 








JCMT 


328 - 366 


0.625 


0.51-0.57 


ACSIS 


4-9 


90-400 


14 


0.56-0.53 


18 


26 


25 



(1) In the 2mm band, each frequency setting is observed twice, i.e. once at each of the two spectral resolutions (0.32 and 1 MHz), while in the 
197-265 GHz band, each frequency setting is observed only one time, at a slightly different spectral resolution (1 or 1.25 MHz), (2 'rms is given in 
1.5 km.s~' bins, {i) crcai is the calibration uncertainty, (4 'N„,; is the number of compared lines for calibration purposes, <5) P ra / is the percentage of 
compared spectra for calibration purposes, (6 'these values refer to the "internal" comparison between VESPA and 1 MHz FB spectra simultaneously 
observed; the other ones refer to "external" comparisons with published or previously obtained observations. 



survey spectra that could be cross-checked and calibrated against 
previously published data for each band. In addition, we cross- 
checked the calibration in the 2 mm band by comparing the data 
obtained with 1 MHz and 320 kHz resolution. Finally, we esti- 
mated the calibration uncertainty due to the receivers by com- 
paring lines in the 1 mm band observed with the two receivers 
A230 and B230. 

To quantify the differences, we obtained Gaussian fits of 
the considered lines and compared their characteristics (inte- 
grated intensity, peak intensity and full width at half maximum, 
FWHM), with the previously published values. Note that the 
comparison was performed in the main beam brightness scale 
(T,„i), based on the T^ fT m b beam efficiency factors given in 
Table[T] This method provides two types of check: i) the average 
uncertainty over each band; and ii) possible specific calibration 
problems on single settings. 

3.2. Calibration uncertainties 

With the method described above to quantify the calibration 
uncertainty in the survey, we obtained the following results. 

FWHM: 

In all the frequency ranges, except 1mm, the agreement between 
the survey and the published FWHM values is within 15-20% 
and no systematic trend is observed. In contrast, in the 1 mm 
band the FWHM of the survey lines appears systematically 
broader by 1 kms -1 than the published values. This is likely 
due to the relatively poorer spectral resolution of our survey in 
this range (0.8 to 1.9kms _1 ) compared to the linewidths (on 
average ^ 4-5 kms -1 ). 

Integrated and peak intensity: 

The comparison of the integrated and peak intensity of the lines 
of the survey with published values yields the same results, 
when the difference due to the spectral resolution described 
above is taken into account. Besides, since the derivation of 
the integrated intensities does not depend on the line shape, we 
choose to quantify the calibration uncertainties by comparing 
the integrated line intensities. 

Derived uncertainties: 

We can derive an estimate of the calibration uncertainty from 
the distribution of the integrated intensity ratios. Fig[T]shows the 
result for the 1mm band. If one excludes the two extremes at 



< 0.5 and > 1 .5 which correspond to spectra with "anomalous 
intensity", the ratios distribution can be fitted by a Gaussian with 
a mean value R, very close to 1 and a standard deviation o~ ratw . 
With the assumption that the relative uncertainties on the pub- 
lished intensities, cr pu i,, and on the survey line intensities, cr ca i, 
are independent variables, the error propagation formula leads to 
the following relationship: 

<r 2 rano=R 2 (<rl aI + CT 2 pllh ) (1) 

Most publications report calibration uncertainties of 15% (or do 
not report any estimate). Except in the 3 mm frequency range, 
where our comparisons suggest that the calibration uncertainty is 
probably somewhat lower than 15%, we obtain consistent uncer- 
tainties either assuming cr pu b = 15% for all the frequency ranges 
or assuming that our observations are representative of average 
observation conditions in each frequency range, i.e. cr pu b = <x cfl/ 
and thus cr^ a[ ^ cr^ fl(( . o /2. Table fllreports the resulting calibration 
uncertainties for each observeaband. 

It can be noticed that at 2 mm, the "external" comparisons 
with published spectra, which include all uncertainty factors, 
lead to higher values than "internal" comparisons between the 
VESPA and 1 MHz FB simultaneous observations, which take 
into account only the contribution of the backends. Similarly, 
at 1 mm, comparison of the line intensities observed simulta- 
neously with the two receivers A230 and B230 shows that the 
receivers' contribution to the calibration uncertainties is of the 
order of 10%, whereas the total calibration uncertainty is 17%. 

Finally, the full list of lines used for the calibration compar- 
ison is reported in the On Line Material, in Table [2] for a com- 
parison between the survey lines and previous observations ob- 
tained with the same telescopes and in Table[3]for a comparison 
between the survey lines observed simultaneously with VESPA 
and the 1 MHz Filter Bank at the IRAM 30 m telescope. 

4. Data release 

The data are made publicly available on the web site http://www- 
laog.obs.ujf-grenoble.fr/heberges/timasss The site provides the 
files with the IRAM-30m and JCMT- 15 m d ata in CLASS for- 
mat (h ttp://www.iram.fr/ IRAMFR/GILDAS ). The intensities are 
in Ta*. Based on the discussion of the previous section, the 
potential user is highly recommended to verify the calibration 
uncertainty of the data that she/he wants to use by looking at 
Tables|2]and|3] We emphasize that we did not apply any "rescal- 
ing" factor to the survey data because the difference may be 
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Fig. 2. The IRAS 16293 spectra in the four bands of the survey. Upper panels: IRAM-30m 3 mm and 2 mm bands. Lower panels: 
IRAM-30 m 1 mm and JCMT-15 m 0.8 mm bands. The middle panels are blow-ups of sample frequency ranges in the four bands 
respectively. These panels include lines identification based on the publicly available spectral databases (see text for details). 
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Fig. 1. Distribution of integrated intensity ratios of this survey's 
1 mm lines compared to published observations obtained with 
the IRAM 30 m telescope. The curve is a Gaussian fit to the 
histogram ruling out the "anomalous" ratios (< 0.5 or > 1.5). 
It can be noticed that the 12 "anomalous" ratios correspond to 
only 7 "anomalous spectra" among the 165 spectra observed in 
the 1 mm range. 



caused by a wrong calibration of the published data rather than 
a wrong calibration of the survey data. Only a very careful sci- 
entific analysis can assess what is the best. It is, therefore, the 
user's responsibility to verify that the data are correctly used, 
based on the information provided on the web site. 



5. Results 

5.1. Overall survey 

Figure|2]shows the full survey in the four bands and the richness 
of the IRAS 16293 line spectrum. On average about 20 lines per 
GHz have been detected with a S/N larger than 3 in the 220 GHz 
frequency range covered by the survey. The line density seems to 
slightly increase with frequency: 17 lines/GHz in the 3 mm band, 
19 and 23 lines/GHz in the 2 mm and 1 mm ranges respectively 
and up to 26 lines/GHz in the 0.9 mm range. 

To quantify more rigorously the lines and species detected in 
the survey, we made Gaussian fits and line identification using 
the CASSIS package (http://cassis.cesr.f^). The spectroscopic 
data come from the CDMS and JPL databases ( jMiiller et al.| 
[2001] [2005] |Pickett et al^[l998l and references quoted on the 
databases to data producers for each species). In a few cases (or- 
tho and para H2CO for instance) a specific database with each 
form separated has been used (see http://cassis.cesr.fr). For the 
D-bearing isotopologues of methanol, only the lines reported by 
( Parise et al.|20 02 2004) are included in this paper. 

Hereinafter we will only consider lines identified according 
to the following criteria: i) lines belonging to species included 
in the JPL and CDMS databases or to the D-bearing isotopo- 
logues of methanol, ii) lines detected with more than 3<x in the 
integrated line intensity, iii) unblended lines and iv) lines with 
upper level energies E up lower than 250 K. This last condition 
only limits the number of methanol lines in this analysis, since 
such lines of other molecules are too weak in any case. When ap- 
plying these criteria, we end up with ~ 1000 lines listed in Table 
|4](Online Material). In the table we report the line identification 
together with the result of the Gaussian fit of each line (see also 

Figure \M shows the line densities, limited to the lines sat- 
isfying the above criteria, in each of the four survey frequency 
ranges, for various signal to noise ratios. In the 3 mm and 2 mm 



3<S/N 




10 < S/N 



2mm 



0.9mm 



3mm 



1mm 
Frequency ranges 



Fig. 3. Distribution of the density (averaged over 10 GHz inter- 
vals) of identified lines (see text) in each of the four frequency 
ranges. The upper panel corresponds to lines with S/N < 3, the 
central panel is restricted to lines with 3 < S/N < 10 and the lower 
panel corresponds to lines with signal to noise ratios > 10. 



10000 
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CASSIS has been developed by CESR-UPS/CNRS 



o 



0.3 1 3 10 30 

Line integrated intensity (K.km/s) 

Fig. 4. Number of lines showing an integrated intensity larger 
than a given threshold. The solid line shows the single index 
power law best fit of the distribution. 



ranges, theses densities are a factor 2 to 3 smaller than the esti- 
mates of the total densities including blended lines. This effect 
is still stronger in the 1 mm and 0.9 mm ranges, where the lack 
of weak lines (3 < S /N < 10) is particularly striking in Fig. 
[3] This is a bias due to our selection criterion of non-blended 
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lines. In fact, these frequency ranges are rich in lines from large 
molecules, which emit many weak lines, so that our unblended 
line criterium filters out a large fraction of these moderate signal 
to noise ratio lines. In contrast, the 2 mm range, which benefits 
from a better spectral resolution compared to the 3 mm and 1 mm 
ranges, suffers less from this selection effect. 

Overall, most of the lines retained in the 1mm and 0.9 mm 
ranges and a large fraction (~ 2/3) of the lines retained in the 
3 mm and 2 mm ranges have a high signal to noise ratio 10). 
The density of lines with such high S/N is relatively constant in 
frequency and equal to ~ 4 - 5/GHz. 

When comparing with the line density quoted at the begin- 
ning of the section, namely about 20 lines per GHz, obtained 
considering lines with S/N>3 but no other filter, clearly the in- 
troduction of the other criteria, unblended lines and, to a lesser 
degree, identified lines and E up < 250 K, severely underestimates 
the line content. 

The line intensity spans more than three orders of magnitude, 
from 1 mK to 24 K. The number of lines showing an integrated 
intensity larger than a given threshold is given in Fig. [4] For inte- 
grated intensity ranging between 1 and 30 K.km/s, the distribu- 
tion roughly follows a power law of slope -0.9. The power law 
breaks down in the high and low end of the distribution. This 
slope is identical to the those observed by |Schilke et al.| ( |2001) >; 
[White et aL] ( |2003| >; |Comito et aL] ( |2005| l in their submillimeter 
surveys of Orion-KL. Similarly with what is observed in these 
surveys, this slope does not provide a good fit for the brightest 
and the weakest lines. 



5.2. Detected lines and species 

Table [5] lists the species detected and identified in the survey, 
with the number of lines, the range of upper level energy of the 
lines and other observational quantities (see below). 

In the frame of a survey analysis restricted to the four 
line criteria mentioned above, 69 different molecules (includ- 
ing ions) have been detected. They correspond to 32 distinct 
chemical species and include 37 rare isotopologues. Out of the 
~ 1000 lines of Table[5] about half belong to only three species: 
CH 3 HCO, HCOOCH3 and CH3OH. 

Most of the 4000 lines detected in our survey belong in 
fact to the species already identified in this source (see Table 
5), among which a few molecules show extremely rich spec- 
tra, with many weak and/or blended lines, not included in the 
present study. Although we anticipate the presence of unidenti- 
fied lines, their identification will require a careful analysis and 
even probably modeling of these spectra. Also, with the E up < 
250K selection criterion used in this work, we do not identify 
any vibrationally excited lines. These objectives are postponed 
to a future article. 

Table [5] also lists the sum of the line flux in each species. 
For species showing rich and complex spectra or with numer- 
ous lines with upper energy levels > 250 K, this value is in fact 
a lower limit due to the numerous weak and blended lines not 
included in the present analysis. 

As noted by the previous study (|Blake et al. 1994 



Dish oeck et al.|[l995) , the millimeter and submillimeter spec 
trum of IRAS 16293 is dominated by simple O-rich species like 
CO, SO, H 2 CO, S0 2 , and CH 3 OH (the three first families alone 
contribute to more than 2/3 of the total flux). In the frequency 
range covered by our survey, the total flux emitted in the CO 
main and isotopic lines is about 800 K.km.s -1 , ~ 30% smaller 
than the flux emitted in the SO, H 2 CO, S0 2 , and CH3OH lines 



together (~ 1100K.km.s~ 1 ). Thus CO is not the major cooling 
agent in this frequency range. In addition, our survey shows the 
presence of thousands of weaker lines from heavier and more 
complex molecules not detected in the previous surveys. The 
blow-ups of Fig.|2]illustrate the situation, with several lines from 
CH3CCH, CH 3 CN, HCOOCH3 etc. To have a rough estimate of 
the contribution of this "grass" of lines to the total cooling of 
the gas (in this range of frequencies), consider a line density of 
about 15 lines per GHz (see above) and a line integrated inten- 
sity of ~ 0.3 K.km.s -1 . This gives approximatively an integrated 
line intensity of 900 K.km.s -1 , namely comparable to the contri- 
bution of CO and its isotopologues. 



5.3. Line parameters 

A Gaussian fit has been performed for each of the lines of Table 
|4] The parameters of the fit (integrated intensity, line width 
FWHM and rest velocity Vlsr) are reported in the same ta- 
ble. The uncertainties on the integrated line intensity take into 
account the spectra noise and the calibration uncertainties re- 
ported in Table [T] Note that we have verified that even when 
the line profile differs from a simple Gaussian (for example in 
self-absorbed lines, or lines with broad wings) the Gaussian fit 
area is very close to the true integrated line intensity and, there- 
fore, reliable. The uncertainties on the line width FWHM and 
rest velocity Vlsr take into account the statistical errors (from 
the fit) and the uncertainty due to the limited spectral resolution, 
which indeed dominates the error. The cases when the line pro- 
files are clearly not Gaussian are marked in Table|4] The average 
<FWHM> and <Vlsr > for each species and isotopologues are 
reported in Table [5] except for those which show severe blend- 
ing due to unresolved or partially resolved hyperfine structure 
(CH3CN, HNCO, NH 2 D, NS). 

In the following, we analyze the information derived from 
the Gaussian FWHM and Vlsr of the lines. 

Figure [5] shows the FWHM versus v L sr for each species. 
In the plots, we have regrouped the isotopologues of the same 
species and, in a few cases, species with a small number of lines. 
Note that in these cases we verified that the species have simi- 
lar FWHM and Vlsr to avoid introducing artificial trends in the 
plot. Figure [5] shows a remarkable and unexpected behavior: the 
FWHM and v L sr distributions are not the same but, on the con- 
trary, they depend on the species. 

Based on the different FWHM and v l $r distributions, it is 
possible to identify four types of "kinematical behaviors": 

1. Typel (first row of Fig. Bji: <v L sr> ~4.0km.s~ 1 and 
<FWHM> ~ 2.5 km.s -1 . The lines show very little dis- 
persion both in terms of rest velocities and in term of 
line widths. Small carbon chains and rings, and "small 
molecules" belong to this group. 

2. Type II (second row): <Vlsr> ~ 3.7 km.s -1 with a very little 
dispersion, FWHM from ~2 km.s -1 to ~8 km.s -1 . All species 
in this group are S- or N- bearing molecules. It can be no- 
ticed that HCO + and C3H2 (in the first row) show a behavior 
intermediate between Type I and Type II. 

3. Type HI (third row): <v LSR >~ 2.5-3.0 kms -1 , <FWHM> < 
4.0 km.s -1 . Four complex organic molecules show this be- 
havior. 

4. Type IV (fourth row): Vlsr and FWHM showing a mixed be- 
havior, with characteristics belonging to the two last groups. 
CH3OH lines have Vlsr and FWHM ranging from 2 to 
9 km.s -1 ; H 2 CO and CH3CCH lines have moderate FWHM 
(4—5 km.s" 1 ) and Vlsr ranging from 2.5 to 4km. s -1 ; lines 
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Table 5. Detected species ordered by decreasing fluxes of the main isotopologue. 



<V LSR > <FWHM> j T mb dv 
(km.s-') (K.km.s" 1 ) 



<V LSR > <FWHM> j T mb dv 
(km.s- 1 ) (K.km.s- 1 ) 



TAG 

(i) 



Species 

(3) 



Nb< 2 » 
lines 



E up (K) 
min-max 



TAG 

(i) 



Species 

(3) 



Nb (2 > 
lines 



E up (K) 
min-max 



28503 
29501 
30502 
29006 
48501 
50501 
30591 
30581 
31501 
31503 
32592 
32582 
32503 
64502 
66002 
32504 
33502 



29507 
30504 
30510 
31506 
31508 
30505 
44501 
46501 
45501 
45502 
27501 
28002 
28509 
28003 
44003 



(aa) CO 


3 


5.5- 


•33.2 


5.93±0.71 


9.47+1.77 


605.8 


(fl 6)13 CO 


3 


5.3- 


•31.7 


3.94±0.36 


4.13+1.16 


116.0 


( flc ) c 18Q 


3 


5.3- 


31.6 


3.72+0.18 


2.90+0.61 


67.6 


<w> C "o 


3 


5.4- 


•32.4 


3.84±0.35 


3.75+0.66 


25.5 




20 


9.2- 


•87.5 


3.98±0.20 


4.18+0.55 


339.5 


( fl /)34 SO 


13 


15.6- 


86.1 


3.81+0.27 


4.87+0.84 


28.0 


^o-HjCO 


10 


6.8- 


■143.3 


3.45+0.58 


4.09+0.58 


163.3 


(<,s) p-H 2 CO 


7 


10.5- 


99.7 


3.67+0.12 


3.92+0.65 


93.6 


^HDCO 


10 


17.6- 


•102.6 


3.39+0.58 


3.99+1.06 


26.7 


w H 2 13 CO 


12 


10.2- 


157.5 


2.96+0.38 


2.86+1.29 


12.7 


(< ^o-D 2 CO 


6 


16.8- 


81.2 


3.51+0.61 


3.80+0.84 


11.7 


( ^p-D 2 CO 


8 


5.3- 


•99.5 


3.38+0.72 


2.76+1.24 


5.5 


(*H 2 C 18 


4 


20.0- 


-63.4 


2.71+0.16 


2.67+1.21 


1.9 


(fl,) so 2 


82 


7.7- 


•248.5 


3.74+0.32 


5.30+1.22 


238.3 


(am) 34 gQ^ 


5 


7.6- 


69.7 


3.33+0.72 


4.20+1.71 


1.6 


•"'"CHsOH 


99 


7.0- 


•233.6 


3.42+0.42 


5.91+1.35 


217.6 


(flo)13 CH 3 OH 


9 


6.8- 


•94.6 


2.93+0.47 


3.92+1.63 


2.8 


^CHjDOH 


12 


4.5- 


-40.8 


2.74+0.64 


5.89+0.91 


6.4 


(flp) CHb 2 OH 


7 


4.2- 


-53.5 


3.70+0.88 


4.76+1.58 


2.6 


(flp) CH 3 OD 


4 


6.0- 


•26.8 


2.55+0.84 


3.96+1.29 


2.3 


<"'i ,) CD 3 OH 


3 


12.2- 


•42.0 


2.43+0.21 


2.03+0.49 


0.5 


(ar) HC Q+ 


3 


4.3- 


•42.8 


3.29+0.16 


2.87+1.39 


144.2 


( flr ) H 13 CO + 


4 


4.2- 


41.6 


3.94+0.46 


2.43+0.44 


31.5 


( " s) DCO + 


3 


10.4- 


51.9 


4.19+0.29 


2.34+0.83 


14.3 


(fl,) HC 18 + 


3 


4.1- 


•24.5 


4.04+0.37 


2.15+0.51 


1.6 


(ar) D 13 CO + 


2 


10.2- 


•20.4 


4.52+0.16 


2.34+1.36 


0.6 


<"" ) HC 17 + 


1 


12.5- 


•12.5 


3.94+0.22 


4.03+0.71 


0.4 


<ar) CS 


4 


7.0- 


•65.8 


3.70+0.12 


3.56+0.43 


126.5 


(a ,) C 34 S 


4 


6.9- 


•64.8 


3.63+0.16 


3.60+0.68 


21.1 


(ov)13 CS 


4 


6.7- 


•46.6 


3.69+0.22 


4.06+0.96 


9.0 


(ov) C 33 S 


4 


7.0- 


•65.3 


3.68+0.66 


5.05+1.67 


7.2 


(a,v) HCN 


3 


4.2- 


•42.5 


3.57+0.78 


7.06+2.12 


97.5 


<M) H 13 CN 


5 


4.1- 


41.4 


3.77+0.26 


4.05+1.85 


17.1 


fe)DCN 


3 


10.4- 


52.1 


3.63+0.46 


4.09+0.55 


11.3 


te) HC I5 N 


4 


4.1- 


41.3 


3.60+0.33 


5.52+1.47 


6.4 


(Afl, CH 3 CHO 


192 


13.8- 


•194.5 


2.79+0.39 


2.68+1.00 


70.0 



41001 
46509 
47505 
47504 
60503 
62505 
61502 
60003 
44505 
27002 
28508 
28005 
51501 
52005 
34502 
35001 
38082 
38092 
39003 
46008 
40502 
41502 
25501 
26501 
43002 
26504 
45506 
42501 
18501 
30008 
56007 
45013 
19002 
46010 
49503 
37003 



(W "CH 3 CN 


>30 


20.4- 


•21.2 






58.6 


{bc) H 2 CS 


26 


9.9- 


■143.4 


3.61+0.35 


4.53+1.34 


51.2 


(M) H 2 13 CS 


4 


22.5- 


97.1 


4.09+0.79 


4.61+2.29 


2.9 


•'""'HDCS 


5 


8.9- 


■30.9 


3.77+0.12 


2.81+0.88 


0.8 


W'OCS 


12 


16.3- 


■134.8 


3.13+0.33 


5.12+0.77 


38.6 


(* s ) OC 34 S 


9 


15.9- 


•157.2 


2.51+0.26 


3.98+0.65 


7.8 


(i s ) 13 CS 


10 


20.9- 


■147.2 


2.58+0.51 


4.49+1.25 


5.6 


(W,) HCOOCH 3 


182 


17.9- 


•145.0 


2.52+0.32 


2.71+0.77 


51.6 


(W) SiO 


6 


6.2- 


75.0 


4.40+0.46 


5.57+0.78 


49.7 


(«HNC 


3 


4.3- 


43.5 


3.80+0.56 


2.81+1.08 


32.1 


(M) DNC 


2 


11.0- 


•22.0 


4.31+0.02 


2.47+1.19 


4.9 


whn 13 c 


4 


4.2- 


41.8 


3.74+0.61 


3.44+1.76 


4.8 


(W) HC 3 N 


17 


19.6- 


•203.0 


3.53+0.44 


4.67+1.97 


35.4 


''""'DCsN 


6 


22.3- 


77.0 


4.29+0.23 


2.61+0.41 


1.0 


(i "°H 2 S 


2 


27.9- 


84.0 


3.55+0.52 


4.90+0.08 


25.2 


(M H ds 


2 


11.7- 


34.7 


3.05+0.55 


2.41+2.09 


0.7 


^o-c-CjHz 


14 


4.1- 


•48.3 


4.10+0.27 


2.73+0.70 


14.5 


^"'p-c-CjHz 


9 


6.4- 


47.5 


4.27+0.26 


2.57+1.29 


4.5 


(fo?) c-C 3 HD 


5 


17.4- 


•25.34 


4.44+0.31 


1.91+0.88 


1.5 


(6r) CH 3 OCH 3 


82 


11.1- 


169.8 


2.53+0.29 


2.27+0.69 


18.8 


(i,) CH 3 CCH 


26 


12.3- 


•176.6 


3.62+0.32 


3.03+0.58 


17.4 


< 6,) CH 2 DCCH 


4 


16.3- 


38.1 


4.34+0.31 


1.93+0.40 


0.2 


(i„) CCH 


13 


4.2- 


25.1 


3.83+0.09 


2.08+0.20 


14.7 


(*v) CCD 


4 


10.4- 


10.3 


4.32+0.52 


2.54+0.86 


0.8 


(6,r) HNCO 


>8 


10.5- 


•143.5 






13.8 


(*-v) CN 


7 


5.4- 


5.4 


4.45+0.10 


1.82+0.16 


8.5 


(6v) HCS+ 


5 


6.1- 


73.7 


3.72+0.20 


3.32+0.70 


5.1 


(i;:) H 2 CCO 


13 


9.7- 


•66.9 


3.00+0.26 


2.58+0.44 


4.8 


( '' fl) NH 2 D 


>3 


16.6- 


•182.8 






4.7 


<<*>NO 


7 


7.2- 


36.1 


4.02+0.12 


2.25+1.42 


3.7 


( '' f) C 2 S 


12 


15.4- 


•53.8 


4.10+0.17 


2.12+0.62 


2.6 


(£,rf) PN 


>3 


13.5- 


63.1 






2.1 


{ce> HDO 


1 


95.2- 


•95.2 


2.87+0.21 


6.13+0.58 


2.1 


(C/) NS 


>6 


9.9- 


17.8 






1.3 


(CS) C 4 H 


6 


20.5- 


•30.2 


4.00+0.08 


1.58+0.28 


0.7 


( <*>c-C 3 H 


3 


4.39- 


10.8 


4.42+0.22 


1.79+0.10 


0.3 



(1) In the CDMS as in the JPL catalog, the species are associated to a five digit number, called TAG; the first two digits correspond to their 
molecular weight in atomic mass units, and the third digit is a 5 in the CDMS catalog, and a in the JPL catalog. Our identification of the ortho and 
para forms of H 2 CO, D 2 CO and c-C 3 H 2 relies on the VASTEL Spectroscopic database (http://www.astro.caltech.edu/~vastel/CHIPPENDALES) 
so that these species are associated with specific TAG. The deuterated isotopologues of methanol, not yet included in the CDMS and JPL 
databases, do not have a TAG. Our identification of their lines is based on the data given in Parise et alT| {2002| and |Parise et al.||2004[ > and 
references therein. 

(2) For most N-bearing species, due to unresolved hyperfine structure, this number is a lower limit corresponding to groups of blended lines. 

(3) For each species, the spectroscopic refe rences given here are the m o st recent cited i n the C DMS and JPL databas es. All of them should b e 
read as reference and references therein: CK Winnewisser et al. Ul997| >; jCazzoli et all (12004 ); jKlapper et all^OOl) ; HKlapper et ~aT1 (|2003|>; 
"iBogey et a L| | l997); "{Klaus et al.|(|l996|; "^MUller et a l. (2000d); "'Bocquet et al. (1999); " Muller et al. (2000c); " Lohilahti & Horneman 
(2004); M Muller et al. (2000b); a Mul ler et al.| ( |2000a|>; fl "|Belov e t al.|^1998>; "IMuller et al.|(f2004); "jXu & Lovas||T997|); ^Parise et aL](|2002|>; 
"•jParise et al.|(|2004|>; "ILattanzi et al.|(|2007t; "Ivan der Tak et al.|d2009t; ISchmid-Burgk et al.U2004|); "1Dore et al.U200U; ^i5m & Yamamoto 



( 2003 ); "" Thorwirth et al. ( 2003 1; c " Cazzoli & Puzzarini ( 2005 1; fll Briinken et al. ( 2004 ); "■ Cazzoli et al. ( 2005 ); 
let all |2009|>; "IClouthier et al.|dl994|>; IBrown et al.H1987>; " I Minowa et al.H1997>; "feolubiatnikov et al ' " 
Maeda et al. ( 2008 1; "|Cho & SaTtoT(|1998|); ^Thorwirth et al. ( 2000a I; " k Brunken et al.|l2006 1; *' Thorwirth et al 



' Kleiner et al. 



(1996); 



'Muller 



200 5>; ^Lovas^ Suenram 
(2000b); ""■Spahnetal. 



1987 



2008 



'" Belov ( 1995 ); Camy-Peyret et al. ( 1985 ); *' Mollaaghababa et al. ( 1993 ); "''Bogey et al. ( 1987 ); "'Neustock et al. ( 1990 1; '"Cazzoli & Puzzarini 
(2008); fa Legue nnec et al.lfliggSllf^Pa d ovani eTa l. (2009l;' 5 fB ogev et al.||l985|>; "ILapino v et al.|(|2007|>; "'Klisch et aTTfl995k "'Margules et al 
(2003 • b: 



•Guarnieri & Huckauf (2003); "Fusina et al. ( 1988); ''Meerts & Dymanus ( 1972); "Lovas et al. 1 1992); "Cazzoli et al. (2006); 
tl983l ; c ^Lovas & Suenram| ( |1982| iT lc"hen et al.U1995| >; c 1Yamamoto & SaitoH1994 ; 



'Johns 



from the rare isotopes of OCS have Vlsr ~ 2.5km.s 1 and 
FWHM<4km.s I . 

It can be noted that the FWHM averaged over each of the 
four frequency bands of the survey is similar, between 2 and 
5 kms _I , with a slight increase in the 1 mm band due to a poorer 



spectral resolution. We have verified that none of the four kine- 
matical behaviors is an artifact due to this instrumental effect. 
Table [6] summarizes the situation. In order to better understand 
the physical meaning of the plots in Fig.|5](and of the four iden- 
tified Types) we have plotted in Fig.|6]the values of the FWHM 
as function of the upper level energy E up of the transition. The 
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species have been grouped as in Fig. [5] it is striking that the dis- 
tinction between the four Types defined by the (FWHM, Vlsr) 
distribution is visible also in this plot. 



Table 6. Distribution of the detected species in four kinematical 
types. 



Type 
(Source) 


v. ™ FWHM 
* LSR r vv lllyl 

(kms 1 ) 


p 

(K) 


FWHM 
r vv jnivi 

behavior 




Type I 
(Envelope) 


~ 4 < 2.5 


0-50 


constant 


CN, NO, 
C 2 S, C 2 H, 
C^H, C 4 H, 
HCO + , C 3 H 2 


Type II 
(Source A) 


-3.7 2-8 


0-250 


increases 


HCN, HC 3 N, 
HNC, 

so, so 2 , CS, 

HCS+, H 2 CS 


Type III 
(Source B) 


2.5-3 < 4 


0-200 


constant 


CITCHO, 
HCOOCH 3 , 
CH 3 OCH,, 
H 2 CCO 


Type IV 
(mixed) 


2.5-4 2-8 


0-250 


increases 


CH,OH, 
ITCO, 
CH3CCH, 
OCS 



First column reports the associated component (see ^6j, second, third 
and fourth columns report the typical velocities, FWHM and upper level 
energy E up ranges of the detected lines. Fifth column describes the be- 
havior of FWHM with increasing E up and last column lists the species 
belonging to each type. 



to the higher sensitivity (~ lOmK versus ~40mK). The average 
line density of the previous unbiased survey of IRAS 16293 was 
7 lines per GHz to compare with 20 lines per GHz for the present 
one. Several species detected in our survey were not detected 
in the previous one: complex organic molecules (HCOOCH3, 
HCOOCH3 and CH3OCH3), carbon chains and rings (C 2 S, C 4 H, 
C-C3H), N-bearing species (NO, PN, NS) and several D-bearing 
molecules. 

Towards hot cores, numerous surveys have been performed. 
They cover the whole range of frequencies reachable from the 
ground, from the 3 mm range observed with IRAM, SEST, 
NRAO or JCMT to the submillimeter windows observable with 
the CSO and JCMT dMacDonald et al.|Jl996; Schilke et at] 



1997; Helmich & van Dishoeck 1997; Nummelin et al. 1998 



[ Thompson & MacDonald|[T999t [Kim et al.|[2000t |Nummelin 
et al.1120001 |Lee et al.||20011 |Scnilke et al.||2001t |White et al 



2003 i lComito et al.|20051|Kim et al.|20061|Belloche et al.|2 007 ; 
Olofsson et al. 2007[ |Tercero e t al. 2010Jl7The line densities usu- 
ally range from 10 to 20 lines per GHz, i.e. comparable to the 
present survey. SgrB2 appears as a noticeable exception, with 
significantly higher line densities, as high as 100 lines per GHz 
in the 3 mm range observed with IRAM ( |Belloche et al.| 2007) 
or 40 lines per GHz in the 1 mm range observed with SEST 
( |Nummelin et al.||T998] |2000| l. Interestingly, the slope of -0.9 



that we observe for the cumulative number of lines versus flux 



threshold is identical to the slopes observed by |Schilke et al. 
d2lM) ; |WhTte et al.| ( |2003] i; |Comito etlu^pOOBT ) in their submil- 
limeter surveys of Orion-KL. 

In conclusion, in terms of molecular content, our survey re- 
veals a richness comparable to that of hot cores and confirms 
that the high abundance of deuterated isotopologues, which are 
easily detected for a number of species, is a distinctive charac- 
teristics of low mass protostars. 



Type I species have lines with E up lower than 50 K. It should 6 - 2 - Kinematical types and associated components 



-up 

be noted that this is not an observational bias: excepted in a few 
cases of very light molecules, such as C2H, the species associ- 
ated to Type I present high energy transitions in the survey fre- 
quency range, but the line intensities decrease abruptly when E up 
becomes larger than 50 K. In contrast, Type II, III and IV species 
present lines with E up up to 200 K but show different behaviors. 
For Type II species, the FWHM increases with E up , whereas for 
Type III and IV the FWHM is constant and does not depend on 
E up . In contrast, analogous plots of v l $r vs E llp show that the 
lines' velocity does not depend on E up in any speci es. A related 
effect has already been observed by Schoier et al. ([2002 ), who 
noted a correlation between the linewidths and the excitation 
temperatures derived by Blake et al. ( 1994| l and |van Dish oeck 
|et al.| ( |1995| l. These properties will be used to give a physical 
meaning to the four Types in ^6] Finally, it can be noticed that 
when detected, the deuterated species show the same behaviour 
as the main isotopomers, except the D-isotopomers of Type III 
species that anyway present too weak lines to be detected in our 
survey. 



6. Discussion 

6. 1 . Comparison with previous surveys 

When compared to the previous survey toward IRAS 16293 
(Bla ke et al.|| 1994 van Dishoeck et al. 1995 1, the present one 
not only enlarges the covered frequency range (~ 200 GHz ver- 
sus ~ 40 GHz) but also the number of detected species, thanks 



As mentioned in the Introduction, IRAS 16293 is formed by a 
pro to-binary system surrounded by an infalling envelope. In ad- 
dition, multiple outflows originate from the system (e.g. [Castets 
|et al. |2001| l. Before attempting to interpret the observations 
of §5.3| and, specifically, the meaning of the four kinematical 
types of Table [6| based on the line rest velocities and widths, 
we review what is known so far about the envelope and the 
proto-binary system. 

a) Envelope 

The envelope extends for 6000-7000 AU in radius, equivalent 
to about 100" in diameter, and it is relatively massive (~ 2M Q ) 
(Crimier et al. 2010| l. At the border of the envelope the dust 
temperature is ~13K and the density is ~ 10 5 cirT 3 . The 
dust temperature reaches 100K at a radius 75-86 AU, where 
the density is (2-3)x 10 8 cirT 3 , creating the region called hot 
corino. The rest velocity of the cold envelo pe has been measured 
in several studies an d it is ~3.9km.s~ 1 (|Mizuno et aL"1 |l990; 
Bottinelli et al.|2004[ |Takakuwa et al.|2007[ >. Molecules probing 
the cold envelope have ~ 2km.s _1 line widths. 

b) Proto-binary system 

Several interferometric studies have been carried out in the past 
to better characterize the nature of the two sources, A (South- 
East) and B (North- West), forming the binary system ([Kuan 



et al. 



(2004); 



Bisschop et a 



Botti nelli et aT| ffiUfy ; |Chandler et aD (|2005 ); 
. (2008 1; see also references in Table |7J. Some 
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v lsr (km.s" 1 ) 



Fig. 5. Plots of the rest velocity Vlsr versus the FWHM, derived from the Gaussian fits of the lines (Table |4|. All the detected 
species and the relevant isotopologues of Table [5] are plotted, except those in which the lines have obviously non Gaussian profiles 
(see text). In particular, the labels HCN* and OCS* mean that the main isotopologues of these species are not included, due to their 
non Gaussian profiles. The 1 sigma errorbars include fit and spectral resolution uncertainties. The vertical lines at Vlsr = 2.7 and 
3.9 km.s 1 correspond to the velocity of the components B and A respectively (see ^6j. The horizontal line at 5 km.s 1 corresponds 
approximatively to the average of the line FWHM range. 



species are associated only or predominantly with one of the two 
sources A and B, and others are observed in both sources. All 
studies agree in measuring broader lines towards A (~ 8 km.s -1 ) 
than towards B (~ 3km.s _1 ). However, the situation about the 
rest velocity of the molecular emission in the two objects is 
more confused: there is some evidence that the two objects have 
different rest velocities (higher in source A than in source B), 
but this difference might also be due to absorption by the en- 
velope according to Chandl er et al.| ( |2005| ) or self-absorption 



et al. | ( |2008] >, both A and B sources show emis sion at velocities 



in optically thick lines from source B according to Botfinelli 
et al. ( 2004 1; Kuan et al. ( 2004} also report velocities somewhat 



higher in source A than in source B, with a dispersion l arger than 
2.5 km.s 1 from one species to another. According to Bisschop 



between 1.5 and 2.5 km.s , Huang et al. ( 2005 ) mention two ve- 
locity components at 1 .5 and 4.5 km.s -1 for source A and show 
emission from B peaking at ~ 2km.s~ 1 . It can be noted that all 
these studies use moderate spectral resolutions (~ 1 km.s -1 ), deal 
with a small number of lines (often only one) for each species 
and, in some cases, suffer from poor signal to noise ratio for the 
weakest lines. 

Table [7] summarizes, for each of the species classified 
according to their kinematical Type in Table [6] the results of 
interferometric observations towards the three components of 
IRAS 16293 (sources A and B and envelope). 
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Type I 



Type II 



Type III 



Type IV 




Fig. 6. Plots of the linewidth, FWHM, versus the energy of the upper level of the transition, E up . The species are grouped as in Fig. 
[5] All the detected species and the relevant isotopologues of Table [5] are plotted, except those in which the lines have obviously non 
Gaussian profiles (see text). In particular, the labels HCN* and OCS* mean that the main isotopologues of these species are not 
included, due to their non Gaussian profiles. The 1 sigma errorbars include fit and spectral resolution uncertainties. The horizontal 
line at 5 km.s~' corresponds approximatively to the average of the line FWHM range. 



c) Interpretation of the kinematical types 

Based on published interferometric observations each species is 
assigned to one or two of the three components (source A, B and 
envelope). For example, CN has the Type I kinematical charac- 
teristics (vlsr ~4km.s -1 and FWHM~ 2km.s~ I ) and has only 
been detected associated with the envelope and was not detected 
in any of the two sources A and B. On the contrary, CH3CHO has 
the characteristics of the Type III behavior (vlsr = 2.5-3 km.s 
and FWHM < 4km.s~ 1 ) and has only been detected in the direc- 
tion of source B. Some species of Table[6](C2H, C3H, C4H, C2S, 
NO, HNC, HCS + and CS) have not been observed with interfer- 
ometers, to the best of our knowledge, so they are not reported 
in Table [7] The correspondence between the four Types of Table 
[6] and the interferometric observations presented in Table [7] sug- 
gests that the species belonging to the same kinematical Type 



are associated with a spatially different source: envelope (Type 
I), source A (Type II), source B (Type III) and a mix of the three 
previous components (Type IV). It should be noted that the dis- 
tribution of molecules in four kinematical types is not an artifact 
or a bias due to the survey pointing, which favors emission from 
source B compared to emission from source A at higher frequen- 
cies. Even excluding all lines observed at IRAM 30 m or JCMT 
telescopes with a HPBW larger than 14", i.e. all lines with fre- 
quencies higher than 200 GHz, the v l $r versus FWHM plots and 
the FWHM versus E up plots show the same behavior. 

1 . Type I corresponds to molecules abundant in the cold enve- 
lope of IRAS 16293. We cannot exclude that these species 
also emit in the densest parts of sources A and B but that 
this emission is strongly absorbed by the gas in the envelope. 
NO, which shows two broad and relatively bright lines with 
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Table 7. Correspondence between kinematical types and spatial 
distributions derived from interferometric observations 



Type 


Source A 


Source B 


Envelope 


Refs. 


Type I 




(C 3 H 2 ) 


CN 


B2004 




HCN 


(HCN) 




T2007 




HCN 






H2005 




HC 3 N 






C2005 




SO 


(SO) 




C2005 


Type II 


SO 






H2005 




so 2 






C2005 




so 2 


(S0 2 ) 




H2005 




H 2 CS 


(H,CS) 




H2005 










K2004 






CH3CHO 




B2008 






CH3OCH3 




C2005 




(CH 3 OCH 3 ) 


CH3OCH3 




H2005 






H 2 CCO 




K2004 


Type III 


(H 2 CCO) 


H 2 CCO 




B2008 






HCOOCH, 




R2006 




HCOOCH3 


(HCOOCH3) 




B2004 




HCOOCH3 


HCOOCH3 




K2004 




HCOOCH3 






C2005 




CH3OH 






C2005 




CH3OH 


CH3OH 




K2004 


TypelV 


H 2 CO 


H 2 CO 




S2004 




H 2 CO 


H 2 CO 




C2005 




OCS 


OCS 




H2005 



The first column reports the kinematical type of the species according 
to the definition given in Table|6] The columns 2, 3 and 4 report where 
the species has been detected: source A, B and envelope respectively. 
The species in parenthesis means that weaker emission has also been 
observed in the relevant component or that they are questionable iden- 
tifications (see text). Last column reports the interferometric observa- 
tions r eferences; B2004:|Bottinelli et al.| < |2004| >; B2008 : BisschopetaL] 
(20081; C2005:|Chandler et al.| < |2005]>; H2005: |Huang et a\.\\2005) -, 
K2004: Ku an et al.N2004^ ; R2006: |Remijan & HollisH2006( i; S2004: 
|Schoier et al.| ( |2004| >; T2007: [ : Makuwa et al.| l |2007| >. 



E up = 36K, bright narrow lines with E up =7K and barely 
detected lines with E up = 19 K (not included in this paper), 
might be an example of such a situation, which deserves a 
detailed excitation study to come. The 343.8 GHz line asso- 
ciated with source B was identified by |Kuan et aL] ( |2004| l 
as the C3H2 231310-2312,11 transition(Ei ow = 790 K) but the 
CH3CHO 182,17 - 172,i6 transition (with E up = 166 K) seems 
a more reasonable identification. 

Type II identifies molecules prevalently associated with 
source A. A few exceptions, for which emission from source 
B is mentioned in the literature, deserve some attention. The 
identification by Huang et al.|(|2005j) of a line to SO2CV2 = 1) 



associated with source B is also questionable given the very 
high upper energy level (Ei ow = 1800 K) of the transition, as 
is the attribution by Huang et al. ( 2005| l of H2CS emission 
to source B, as no map is shown for this species. Finally, SO 
emission is rather extended and associated with both sources 
A and B ( Chandler et al.|2005) . It is possible that part of the 
lines broadening is due to wings associated with the outflow, 
whereas the bulk of the emission is associated with the cold 
envelope (which would make SO rather a species of Type I). 
In addition, it should be noted that, without detailed excita- 
tion models which will be presented in subsequent studies, it 



is not possible to exclude that, for some Type II species, the 
cold envelope contributes to the emission observed at low 
E up and low FWHM. 

3. Type III species have lines prevalently emitted by source B. 
HCOOCH3, which has been observed both in source A and 
source B appears as a notable exception. Excitation mod- 
eling of this species will be discussed in detail in a subse- 
quent paper. Several qualitative arguments may already con- 
tribute to explain why no lines with Type II characteristics 
are present in this study: (i) as the lines emitted by the source 
A are much broader than the lines emitted by the source 
B, those lines present a higher percentage of blending with 
nearby lines which excludes them from this study; (ii) in our 
large beam observations, emission from sources A and B are 
superposed; the narrow lines emitted from the source B ap- 
pear "on the top of broad lines emitted by the source A and 
are, thus, more easily detected; (iii) this second effect is en- 
hanced by the pointing towards the source B. 

4. Type IV species are a mixed bag. The species are present in 
both source A and B and, sometimes, even in the envelope. 
Depending on the intensity contribution, the lines can have 
a low or high v l $r and FWHM. Therefore, Type IV species 
are not associated with any specific component. 

Using the kinematic properties of the lines, we were able to 
propose general trends as whether the emission of the species 
originates from source A or source B although both sources are 
included in our single dish observations. Despite the fact that the 
kinematic differences are small, for example a difference in the 
Vlsr of not more than 1 .5 km.s -1 , we were able to draw a general 
picture thanks to the large number of detected lines. 

6.3. The dynamics of sources A and B 

Having identified where the emission from the various species 
originates helps to clarify the nature of sources A and B. As 
noted in several previous works, A and B have apparently a dif- 
ferent chemical composition: source B is brighter in complex or- 
ganic molecules, while source A is brighter in simpler S- and re- 
hearing molecules (see Table |7J. Furthermore, it was also noted 
that the FWHM of the lines arising in source A are broader than 
in source B. 

Our analysis establishes two additional properties previously 
unrecognized (see discussion in 6.1): i) the two sources have 
different velocities (vlsr) and ii) in the source A, the FWHM 
of the lines increases with increasing upper level energy (Fig.|6]l 
whereas it remains constant in the source B. 

Concerning the linewidths behavior, a related effect, the in- 
crease of FWHM with the excitation temperature, had already 
been observed by Schoier et al. (2002), but this represents a less 
direct probe of the gas kinematics, as it relied on modeling as- 
sumptions to derive the excitation temperatures. The increase of 
the lines FWHM with increasing excitation has been interpreted 
in two ways, either this is due to the effect of the infalling gas on 
the accreting protostellar object ( jCeccareflTe t al. (2003)), or the 
increased width is due to shocks caused by jet/wind interaction 
with the inner dense envelope (Schoi er" et al.| p002), J0rgensen 
|et al.| ( |2002[ )). The two interpretations predict different molec- 
ular emission distributions. However, the existing maps do not 
allow to distinguish these interpretations. In addition, to fully 
disentangle these two interpretations requires a detailed mod- 
eling of the source that is beyond the scope of this paper. At 
this stage, it can simply be said that the interpretation of the in- 
falling gas ( |Ceccarelli et al. 2000b) would lead to reasonable 
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estimates of source A and source B central masses: for optically 
thin emission, the free-fall velocity may be estimated from the 
linewidth, assuming that the FWHM includes quadratically a tur- 
bulent width 6th and the free-fall broadening; its analytic expres- 
sion yields to the following expression of the core masse, where 
the widths are expressed in km.s" 1 and the radius r in AU: 



M(Mq) = 1.4 x IQ- 4 (FWHM 2 - 5 2 th )r 



(2) 



Assuming that the turbulent width in the A and B cores exter- 
nal layers is close to that of the envelope static gas, we derive 
6 t h ~ 2 km.s" 1 from the Type I species linewidths. According to 
the interferometric observations of Type II and Type III species 
reported in Table [7] both sources show a radius of ~ 1.5" (i.e 
180AU at a distance of 120 pc). For core A, typical FHWM 
larger than ~ 6 km.s" 1 lead to a central mass of at least 0.8 M Q , 
whereas for core B, where FHWM are smaller than ~ 3 km.s" 1 , 
the central mass cannot exceed 0. 1M . This assumes that the 
linewidths are due to collapsing envelopes rather than rotating 
disks. If they were, on the contrary, due to rotating disks, one 
would have to take into account the (unknown) inclination an- 
gle. The available interferometer observations, which barely spa- 
tially resolve the core A only (Bottinelli et al. 2004), do not allow 
to distinguish between these two possibilities. In the following, 
we will restrict the discussion to the infalling envelopes, for sim- 
plicity, but our conclusions substantially apply also to the disk 
case. 

The difference between the Vlsr of both sources can be inter- 
preted in the following way: the source B rotates around its com- 
panion, the source A; for a distance of 480 AU between cores A 
and B, and a core A mass of 0.8 M Q , the rotation velocity of 
core B is 1.2km.s _1 . The observed difference between A and B 
velocities (vlsr of 3.9 and 2.7 km.s" 1 respectively) is thus per- 
fectly consistent with this picture. Bottinelli et al. ( 2004[ ) already 
discussed this possibility but did not have enough evidence for 
preferring it with respect to another hypothesis where the lines 
kinematics properties reflect opacity effects. The large number 
of lines and species observed in our survey allows to favor the 
dynamical issue. 



7. Conclusions 

We have presented the unbiased spectral survey of the 3, 2, 1 and 
0.9 mm bands accessible from ground towards the Class source 
IRAS 16293. For that, we used the IRAM-30 m and JCMT-15 m 
telescopes, during about 300 hours of observations. This is the 
most sensitive survey ever published in these bands towards a 
solar type protostar. 

The data have been released for public use in two CLASS 
files, which can be retrieved on the web site http://www- 
laog.obs.ujf-grenoble.fr/heberges/timasss/ The site also contains 
two accompanying files (reported in the On Line Material) pro- 
viding informations about the calibration of the single receiver 
settings, obtained by comparing the survey lines with previously 
obtained observations. 

The line density, ~ 20 lines/GHz, appears as high as in com- 
parable surveys obtained towards high mass protostars (with the 
exception of SgrB2). More than one thousand unblended lines 
with S/N> 3 and upper energy levels lower than 250 K have been 
identified through the comparison with the JPL and CDMS cata- 
logs. They correspond to 32 chemically distinct species, showing 
a chemical richness comparable to that of hot cores. The iden- 
tification of 37 additional rare isotopologues and, specifically, 
numerous D-bearing molecules confirm that IRAS 16293 has a 



remarkably high abundance of deuterated species. The 3 mm- 
0.9 mm spectra are dominated by relatively simple molecules 
(CO, SO, H 2 CO, S0 2 and CH 3 OH). However, the numerous 
weaker lines emitted by larger molecules account for at least as 
much as the CO integrated line intensities. 

The analysis of the profiles of this large set of identified 
lines, and specifically the central velocities and widths, gives 
clues to disentangle where each emission dominantly originates 
from: cold envelope (narrow lines at Vlsr — 3.9 km.s" 1 ), source 
A (broader lines at Vlsr -3.9 km.s" 1 ) and source B (narrow 
lines at Vlsr - 2.7 km.s" 1 . Furthermore, in source A, the line 
widths increase with the upper energy level of the transition. 
If this behavior is interpreted as due to gas infalling towards a 
central object, the core A mass is ~ 1 M Q and the lower line 
widths observed towards source B set an upper limit to the mass 
of this source, < 0.1 M . The observed difference in the Vlsr, 
~ 1.2 km.s" 1 , is consistent with the source B rotating around 
the more massive source A. From a chemical point of view, the 
source B shows predominant emission from O-bearing complex 
molecules whereas N- and S- bearing molecules are strong emit- 
ters in the source A. In order to derive reliable estimates of the 
corresponding chemical abundances, it is necessary to carry out 
careful radiative transfer modeling. This is postponed to future 
articles. 
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Table 2. Comparison of the spectral survey lines with previous observations. 



Setting 


Species 


Transition 


Freq. 


J' 


T A dv 


Ref. 


number 






(MHz) 


(K.kms" 1 ) 












survey 


prev. obs. 




2 


HDO 


lin - In 
± 10 ± 1 1 


80578 


0.507 


0.4 


(6) 


11* 


CH 3 CCH 


5 7 - 4 2 


85450 


0.188 


0.18 


(3) 


11* 


CH 3 CCH 


5. - 4i 

•J 1 ^1 


85455 


0.411 


0.5 


(3) 


11* 


CH 3 CCH 


5n - 4n 
^o 


85457 


0.524 


0.54 


(3) 


14* 


H' 3 CO + 


1 - 


86754 


4.351 


3.97 


(12) 


14* 


H 13 CO + 


1 - 


86754 


4.351 


4.56 


(13) 


14* 


SiO 


2- 1 


86846 


1.910 


1.9 


(4) 


14* 


SiO 


2 - 1 


86846 


1.910 


2.04 


(12) 


14* 


SiO 


2 - 1 


86846 


1.910 


2.23 


(14) 


14* 


SiO 


2 - 1 


86846 


1.910 


2.16 


(15) 


19 


HCO + 


1 - 


89188 


15.195 


13.93 


(13) 


19 


CH 2 DOH 


■MJz x 01 w l 


89251 


0.182 


0.19 


(10) 


19 


CH 2 DOH 


^02 1 1 ^0 


89407 


0.262 


0.31 


(10) 


21 


CH3OCHO-E 


cSno — / r\n 
'-'Uo '0/ 


90227 


0.164 


0.15 


(3) 


21 


CH3OCHO-A 


8ng — 7 07 


90229 


0.144 


0.11 


(3) 


36 


34 SO 


3 7 - 2, 


97715 


0.791 


1 


(4) 


38 


CH3OCHO-A 




98682 


0.248 


0.22 


(3) 


38 


CH3OCHO-E 


(S A — / A A 
i -'4j ' 44 


98711 


0.118 


0.14 


(3) 


60* 


c 18 o 


1 - 


109782 


9.900 


9.02 


(16) 


60* 


c 18 o 


1 - 


109782 


9.900 


7.35 


(17) 


60* 


c 18 o 


1 - 


109782 


9.900 


8.7 


(13) 


61* 


CH 3 CN 


"40 --'40 


110349 


0.164 


0.48 


(3) 


61* 


CH 3 CN 


hiri — Jin 


110364 


0.189 


0.88 


(3) 


61* 


CH 3 CN 


r)">A — JlA 

U ZU — VO 


110375 


0.304 


0.77 


(3) 


61* 


CH 3 CN 


Din — Jin / Don — Jnn 


110382 


0.695 


2.02 


(3) 


62 


D 2 CO 




110838 


0.515 


0.64 


(15) 


65 


c 17 o 


1 - 


112359 


2.855 


3.07 


(16) 


65 


c 17 o 


1 - 


112359 


2.855 


2.59 


(17) 


74 


SiO 


3-2 


130268 


1.361 


4.86 


(15) 


74 


SiO 


3-2 


130268 


1.361 


4.61 


(14) 


81 


CH3CHO-E 


' / "Oo 


133830 


0.383 


0.41 


(3) 


81 


CH3CHO-A 


Ion — \ \f\f. 


133854 


0.597 


0.57 


(3) 


81 


CH 2 DOH 


ipii — z.m 1 
' j Uj ^XjL 1 


133872 


0.739 


0.6 


(10) 


81 


CH 2 DOH 


Tot — Z^o 1 O i 


133881 


0.465 


0.34 


(10) 


82 


HDCO 




134284 


1.466 


1.8 


(9) 


82 


HDCO 


2ii - li° 


134284 


1.466 


1.55 


(15) 


85 


so 2 


5 15 - 4 4 


135696 


3.115 


4.6 


(4) 


85 


so 2 


5i5 ^ 4 04 


135696 


3.115 


3.41 


(12) 


88 


H 2 3 CO 




137450 


0.524 


0.58 


(9) 


88 


H' 3 CO 
Hj 3 CO 
CO 


2 - 1 


137450 


0.524 


0.58 


(15) 


88 


2 ' 2 _ j" 


137450 


0.524 


0.76 


(15) 


88 


9 12 _ , n 

^-■12 ^11 


137450 


0.524 


0.76 


(15) 


95 


H 2 CO 


2 1 2 In 


140839 


19.242 


16.6 


C) 


95 


H 2 CO 


^12 1 1 1 


140839 


19.242 


15.57 


(15) 


100 


C 2 H 5 CN 


Inc n — ] jcii / lf)cii — 1 V in 
av O,1Z 1 AV o,ll A ~O,10 


143406 


0.256 


0.55 


(3) 


101 


CH 3 OCH 3 


^3 4 2 — ^2 5 2 


143600 


0.107 


0.42 


(3) 


101 


CH 3 OCH 3 




143603 


0.414 


0.31 


(3) 


101 


CH 3 OCH 3 


1 1 A Cl — / O ^ A 


143606 


0.189 


0.28 


(3) 


146 


CHD 2 OH 


4n — 3n en 


166435 


0.500 


0.48 


(10) 


151 


H 2 S 




168762 


17.511 


17.5 


(4) 


282* 


H' 3 CO 


3 10 - 2°! 


206132 


1.752 


1.6 


(9) 


282* 


CO 


3° - 2* 2 


206132 


1.752 


1.97 


(15) 


282* 


CO 


3° - 2* 2 


206132 


1.752 


2.19 


(15) 


282* 


H 2 3 CO 


3 _ 2 12 


206132 


1.752 


2.48 


(15) 


282* 


SO 


4 5 - 3 4 


206176 


21.367 


36 


(15) 


285* 


CH 2 DOH 




207780 


0.690 


1.01 


(10) 


292* 


H 2 CO 


3l3 - 2i2 


211211 


31.416 


43.8 


(9) 


292* 


H 2 CO 


3l3 - 2i 2 


211211 


31.416 


33 


(15) 


292* 


H 2 CO 


3l3 - 2i 2 


211211 


31.416 


26.4 


(16) 


296 


H 2 3 CO 


3 22 - 2 2 i 


213037 


0.711 


0.45 


(15) 


296 


H 2 3 CO 


3 2 i - 2 2 o 


213294 


0.593 


0.46 


(15) 


303 


H 2 S 


2 2 o _ 2n 


216710 


7.436 


3 


(4) 


304 


SiO 


5-4 


217104 


4.524 


4.83 


(14) 



Table 2. Continued. 
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Setting 


Species 


Transition 


Freq. 




T.dv 




Ref. 


number 






(MHz) 


(K.kms- 1 ) 












survey 


prev 


obs. 




312 


D 2 CO 


4l4-3i 3 


221191 


1.274 




1.6 


TO 


316 


CH 2 DOH 


5 2 3 - 4 M e. 


223071 


0.972 




1.13 


(Id) 


316 


CH 2 DOH 


533 - 4 32 Oi 


223153 


1.204 




0.58 


(10) 


316 


CH,DOH 


5 2 3 - 4 22 ei 


223315 


0.717 




0.38 


(10) 


321 


HDO 


3]2 — 2 2 i 


225897 


2.089 




1.7 


(6) 


322 


CH3OCH3 


1413 


226346 


1.685 




0.89 


(3) 


323 


CH3CHO-E 


12o,12 ~~ 1 lo.ii 


226551 


0.809 




0.49 


(3) 


323 


CH3OCHO-E 


20 2 ,i8 — 19 2 ,i8 


226713 


1.540 




1.19 


(3) 


323 


CH3OCHO-A 


20 2> i9 — 19 2 ,i8 


226718 


0.707 




0.49 


(3) 


323 


CH3OCHO-A 


20,, 19 - 19 U8 


226778 


1.324 




0.84 


(3) 


332 


D 2 CO 


4o4 — 3fl3 


231410 


3.334 




2.43 


(15) 


332 


D 2 CO 


4o4 — 3o3 


231410 


3.334 




2.01 


(16) 


332 


D 2 CO 


4o4 — 3()3 


231410 


3.334 




1.92 


(17) 


332 


D 2 CO 


4o4 _ 3o3 


231410 


3.334 




2.69 


(13) 


332 


D 2 CO 


4o4 — 3o3 


231410 


3.340 




3 


(9) 


332 


CH3CHO-A 


1248 — H47 


231457 


0.520 




0.51 


(3) 


342 


D 2 CO 


4 22 - 3 2 i 


236102 


1.718 




1.86 


(8) 


342 


D 2 CO 


4 22 - 3 2] 


236102 


1.718 




1.9 


(9) 


342 


D 2 CO 


4 22 - 3 21 


236102 


1.710 




1.9 


(15) 


353 


HDO 


2n — 2i2 


241561 


2.268 




2 


(6) 


363* 


HDCO 


4 14 - 3 13 


246925 


2.067 




8.43 


(15) 


363* 


HDCO 


4 14 - 3,3 


246925 


2.067 




7.6 


(9) 


401 


HDO 


2 2 o — 3,3 


266161 


0.214 




0.21 


(6) 


434 


c 18 o 


3-2 


329335 


55.768 




33.6 


(11) 


436 


13 CH 3 OH 


7,7 _ 6,6 


330194 


0.458 




0.51 


(5) 


436 


13 CH 3 OH 


734 - 633 


330408 


1.123 




0.91 


(5) 


437 


13 CO 


3-2 


330588 


55.854 




67.2 


(11) 


446 


HDCO 


5,4-4,3 


335096 


4.550 




3.95 


(2) 


446 


HDCO 


5,4 - 4,3 


335097 


4.550 




4.5 


(9) 


450 


c 17 o 


3-2 


337061 


12.481 




1 1 


(1) 


450 


CH3OH 


3 3 - 4 2 E + 


337136 


0.826 




1.05 


(2) 


452 


H 2 CS 


21,, 10 — 9,,9 


338081 


1.723 




2.4 


(1) 


455 


34 SO 


89 — 7g 


339858 


3.863 




2.98 


(1) 


461 


D 2 CO 


606 — 5fl5 


342522 


1.875 




1.9 


(9) 


461 


CS 


7-6 


342883 


38.535 




51.4 


(II) 


462 


H^ 3 CO 


5,5 - 4,4 


343325 


3.181 




1.37 


(9) 


464 


SO 


8g — 77 


344311 


18.647 




18.1 


(1) 


470 


SiO 


8-7 


347331 


8.546 




5.79 


(1) 


475 


D 2 CO 


6 2 5 — 5 2 4 


349631 


1.506 




1 


(9) 


479 


H 2 CO 


5,5 - 4,4 


351768 


35.945 




38.3 


(9) 


483 


H^ 3 CO 


5 05 — 4()4 


353812 


0.869 




0.6 


(9) 


484 


HCN 


4-3 


354505 


41.068 




63.4 


II) 


485 


H^ 3 CO 


5 2 4 _ 4 23 


354899 


0.557 




0.3 


(9) 


490 


so 2 


154,12 - 153.13 


357241 


3.603 




3.29 


(1) 


490 


so 2 


H48 - H39 


357388 


4.201 




2.25 


(1) 


496 


DCO + 


5-4 


360170 


3.135 




4.02 


(2) 


501 


HNC 


4-3 


362630 


8.264 




11.9 


(11) 


501 


H 2 CO 


5o5 - 4o4 


362736 


18.394 




28.9 


(9) 


503 


H 2 CO 


5 24 - 4 23 


363946 


11.728 




11.9 


(9) 



All observations are obtained with the same telescopes: IRAM 30 m from 80 to 280 GHz, JCMT from 328 to 366 GHz. The coordinates of the 
observed position are cc(2000.0) = 16 ,! 32 m 22.6™, (5(2000.0) = -24° 28' 33"; except for the settings numbers tagged wi t h an asterisk that correspond 



to a slightly different p osition: ffi(2000 0) = 16 * 32^ 22.7'", £(2000.0 ) = -24 
fcazaux et a H d2003|>; fWakelam e t a I1(|2004t; frarise et al .]l |2004| > 



22' 13". Ref erences 
Parise et al.|(|2005a) 



Blake et 



tagged with < 
ai1 jl994t; 



et al.| ( |2000| ir" jParise e't al.H2002) ;" r |Sch6ier et al.| ( |2002^ ,"unpublished IRAM 30 m observations: ' 
""January 1999, '"March 2000T T7 August 2000, ,5 March 1998. 
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Table 3. Comparison of the strong spectral survey lines simultaneously observed with different backends. 
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The two backends used in the 129-175 GHz range were the VESPA autocorrelator (0.32 MHz resolution) and the 1 MHz resolution Filter Banks 
(FB). 
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Caux et al.: IRAS 16293-2422 spectral survey, Online Material p 11 
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